GL-TR -89-0161 


AD-A216  135 


Seismic  Pulse  Broadening  Associated  With 
Fracture  Damage  Caused  by  Explosions  in 
Crystalline  Rock 


Charles  G.  Sammis 


University  of  Southern  California 
Department  of  Geological  Sciences 
Los  Angeles,  CA  90089 


1989  August  30 


Final  Report 

January  1986  -  January  1989 


Approved  for  public  release;  distribution  unlimited 


OTIC 


ELECTE 
DEC  0  7  1989 


% 


GEOPHYSICS  LABORATORY 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 

HANSCOM  AIR  FORCE,  BASE,  MASSACHUSETTS  01731-5000 


la.  REPORT  SECURITY  CLASSIFICATION 
Unclassified 


la.  SECURITY  CLASSIFICATION  AUTHORITY 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


jO 

rvrlT? 

0249  AM.  0704-0  JM 


itt.  o<  Classification  /  downgrading  schedule 


canizatiom  report  numtcfi(S) 


rmi'V'rwn 


U.  NAME  OF  PERFORMING  ORGANIZATION 


University  of  Southern  Calif 


lie.  AOORESS  {City.  SUM,  snd  ZJPCodtJ 


T 


OFFICE  SYM9CL 

(if 


Uepartnumt  of  Geological  Sciences 
Los  Angeles,  CA  90089-0740 


St.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 
DARPA 


Sc  AOORESS  (City,  5ur*.  *nd  Ilf  Coda) 
1400  Wilson  Blvd. 
Arlington,  VA  22209-2308 


•b.  OFFICE.  JYMSOl 
(tt  Mpptkabht) 

NMRO 


1 .  DISTRIBUTION /AVAILA4IUTY  OF  RE 
Approved  for  public  relenso; 
distribution  unlimited. 


S.  MONITORING  ORGANIZATION  REPORT  NUM9(R(S) 

Go-TR-89-0161 


7a.  NAME  OF  MONITORING  ORGANIZATION 
Geophysics  Laboratory 


7b.  AOORESS  (CRy,  Slat*.  tndZIfCodf) 

Hanscom  AFB 

Massachusetts  01731-5000 


9,  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMRER 

F19628-86-K-0003 


10.  SOURCE  OF  FUNOING  NUMIERS 


PROJECT 

TASK 

NO. 

NO 

6A10 

DA 

II.  TITLE  (Mctodb  Sfwtity  OkstiUttion) 

Seismic  Pulse  Broadening  Associated  With  Fracture  Damage  Caused  by  Explosions  in  CrystalliT 


KfiTw  ■ 


12.  PERSONAL  AUTHOR^) 

Charles  G.  Snmmis 


11a.  TYPE  OF  REPORT 

FINAL  REPORT 


1t\  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVERED 

F«OM_L&& _ TO  1/89 


14.  DATE  OF  REPORT  (Ytor.  Month,  Doy)  IIS.  PAGE  COUNT 
1989  August  30  I  98 


COSAT1  CODES 


SUI-GROUP 


I'l.  SUIJECT  TERMS  (Conuwt  on  rov,  to  If  notoatty  and  kMnttfy  by  Mode  numbot) 
damage  mechanics,  strain,  pulse-broadening 


19.  ABSTRACT  (Continue  on  niton*  if  notouory  and  (donbty  by  Modr  nombor) 


Recent  analyse*  of  the  initiation,  propagation,  and  linkage  of  fractures  from  pre-existing  cracks  have  been  used  to 
(emulate  a  damage  mechanics  for  brittle  solids  loaded  in  compression.  The  model  describes  the  accumulation  of  damage  with 
strvin.  It  has  been  used  to  construct  damage  initiation  surfaces  and  failure  surfaces  which  have  been  tested  using  previously 
published  experimental  data.  It  is  shown  that  a  damage  mechanics  rheology  can  produce  seismic  pulse-broadening  observed  in 
underground  tests  in  granitic  rock.  It  is  proposed  that  die  basic  scaling  relation  between  the  size  of  the  pre-existing  cracks  and 
the  stress  requited  to  initiate  damage  "nay  explain  why  models  based  or.  laboratory  data  and  large  scale  laboratory  simulations 
have  failed  to  reproduce  the  pulsc-broadcning  observed  in  the  field.  „ _ ^  /  '  ,y/ 


20.  DISTRIBUTION /AVAILAIILITY  OF  ABSTRACT  121.  AISTRACf  SECURITY  CUSSIFICATIO 

□  UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT.  □  PYIC  USERS  I  Unclassified 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL  122b.  TELEPHONE  (Indudo  Art*  Co do 

James  Lewkowicz  I  (6171  377-3222. 


VI 


OD  Form  1473,  JUN  86 


frtvious  tditlons  trt  obsohto. 


RITY  CLASSIFICATION 
Unclassified 


SEISMIC  PULSE  BROADENING  ASSOCIATED  WITH  FRACTURE  DAMAGE 
CAUSED  BY  EXPLOSIONS  IN  CRYSTALLINE  ROCK 


INTRODUCTION 


Charles  G.  Sammis 
Department  of  Geological  Sciences 
University  of  Southern  California 
Los  Angeles,  CA.  90089-0740 


Aoevsslon  For 


RIIS  GRAfcl 
MIC  TAB 
Unannounced 
Justification. 


t 


□ 


By - 

Distribution/ 


Availability  Codes 
Avail  and/or 
Special 


For  the  three  U.S.  underground  tests  in  granite,  the  particle 
velocity  was  measured  in  the  free  field  as  a  function  of  the  distance 
from  the  shot  point.  In  each  case,  numerical  simulations  using 
continuum  mechanics  with  strength  values  measured  in  the 
laboratory  have  failed  to  predict  the  velocity  pulse  width  to  within  a 
factor  of  about  5  (Rimer  ct  al.,  1987).  Models  which  predict  the 
observed  pulse-width  all  require  a  large  reduction  in  the  post-shock 
strength  to  values  well  below  those  measured  in  the  laboratory  - 
even  for  heavily  pre-fractured  samples. 

The  physical  mechanism  responsible  for  the  post-shock 
weakening  has  not  been  clear.  Cherry  and  Rimer(1982)  invoked  an 
"effective  stress"  mechanism  in  which  the  free  crack  and  pore  space 
is  eliminated  during  the  initial  compression.  Fluids  in  the  remaining 
saturated  cracks  are  compressed  to  high  pressures  and  support  part 
of  the  confining  stress,  thereby  reducing  the  failure  stress  behind  the 
shock-front.  However,  Scholz  (1982)  has  argued  that  dilatancy 
hardening  is  expected  to  override  any  weakening  produced  by  the 
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effective  stress  effect. 

Russian  scientists  (Zamyshlyayev  et  al.,  1980)  have  also  shown 
that  post-shock  weakening  can  produce  the  pulse-broadening 
observed  in  the  Hard  Hat  explosion,  but  their  proposed  mechanism  is 
quite  different.  They  attribute  weakening  to  the  time  dependent 
growth  of  cracks  behind  the  shock  front.  They  also  assert  that  shear 
deformation  behind  the  shock  front  produces  a  rotation  of  the  blocks 
of  rock  between  fractures  and  results  in  dilatancy.  Based  on  their 
modeling,  they  conclude  that  post-shock  weakening  and  dilatancy 
contribute  equally  to  the  observed  pulse  shape.  The  proposed 
weakening  is  shown  in  Figure  1  (from  their  paper). 

Our  explanation  of  weakening  is,  in  principal,  the  same  as  that 
proposed  by  Zamyshlyayev  et  al.  (1980),  except  that  our  model  is 
based  on  the  micromechanics  of  crack  growth  under  compressive 
loading.  As  illustrated  in  Figure  2,  we  model  the  activation  of 
preexisting  cracks  and  joints  (of  length  a)  by  the  explosive  pressure 
pulse.  By  formulating  the  rheology  in  terms  of  the  basic  physics  of 
crack  nucleation,  growth,  and  interaction,  the  effects  of  variables 
such  as  porosity,  water  saturation,  and  the  density  and  orientation  of 
preexisting  cracks  and  joints  can  be  included  in  the  model.  This 
type  of  formulation  is  known  as  a  "damage  mechanics".  The  state  of 
crack  growth  at  any  time  is  represented  as  a  damage  vector.  As 
stress  increases,  damage  accumulates,  the  elastic  stiffness  decreases, 
and  the  stress-strain  relation  becomes  non-linear.  At  large  values  of 
the  damage,  the  growing  fractures  interact,  and  damage 
accumulation  is  accelerated.  A  peak  stress  is  achieved  after  which 
the  strength  decreases  with  additional  strain. 
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In  the  following  we  review  our  damage  mechanics  model 
(Ashby  and  Sammis,  1989;  appended)  and  show  that  it  gives  a  good 
description  of  the  failure  surfaces  of  a  wide  variety  of  rocks  with  a 
very  small  number  of  adjustable  parameters.  We  then  use  the 
damage  mechanics  to  interpret  a  series  of  laboratory  simulations  in 
granite.  Finally,  we  use  a  simple  one-dimensional  calculation  to  show 
how  a  damage  mechanics  rheology  can  produce  significant  pulse 
broadening. 

DAMAGE  MECHANICS 

We  now  briefly  review  our  damage  mechanics  formulation  by 
outlining  the  key  results.  The  details  are  given  in  the  appended 
paper  (Ashby  and  Sammis,  1989)  which  is  in  press  in  PAGEOPH. 

Crack  initiation 

Most  brittle  solids  contain  inhomogeneities:  small  holes  or 
cracks,  particles  which  are  poorly  bonded,  or  phases  which  have 
different  moduli  or  strengths  fiom  those  of  the  matrix.  Any  one  of 
these  can  act  as  nuclei  for  new  cracks  when  the  solid  is  loaded. 

The  range  of  possible  nuclei  is  wide,  but  the  spectrum  of  their 
characteristics  is  probably  bracketed  by  two  extremes:  the  spherical 
hole  and  the  sharp  inclined  crack.  Both  have  been  studied 
experimentally  and  both  have  been  modeled,  the  first  by  Sammis 
and  Ashby  (1986)  and  the  second  by  Nemat-Nasser  and  Horii  (1982) 
and  Ashby  and  Hallam  (1986).  In  both  cases,  the  criterion  for  crack 
initiation,  under  axisymmetric  loading  has  the  form 
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a|  =cl  ^3  ‘  a0 


(1) 


where  cj  and  gq  are  material  properties,  cn  is  the  axial  stress,  and 
02=03  the  radial  stress  (both  positive  when  tensile,  negative  when 
compressive). 

In  this  paper  we  consider  the  growth  ol  crack-damage  from 
initial,  inclined  cracks  as  in  Figure  2.  For  this  case  (Nemat-Nasscr 
and  Horii,  1982;  Ashby  and  Hallam,  1986)  cracks  initiate  when 


where  p  is  the  coefficient  of  friction  acting  across  thi.  crack  faces,  Kic 
is  the  fracture  toughness  of  the  material  through  which  the  new 
crack  propagates,  and  2a  is  the  length  of  the  original  inclined  crack. 
Rocks,  typically,  show  a  coefficient  of  friction  of  about  0.6,  in  which 
case  Cl  =  3.1  and  o0=3.1  K^/ViTa  .  Crack  initiation  from  holes 

(Sammis  and  Ashby,  1986)  gives  similar  values. 

Equation  (2)  was  shown  by  Ashby  and  Sammis  (1989)  to  give 
an  adequate  description  of  triaxial  data  for  a  wide  range  of  rock- 
types.  It  is  used  to  describe  the  initiation  of  damage  in  the  diagrams 
shown  later.  In  each  case,  experimental  data  are  fitted  to  eqn.  (2)  to 
give  p  and  a  (using  published  data  for  Kic).  Results  of  this  analysis 
are  summarized  in  Table  1.  In  the  computations,  it  is  convenient  to 
normalize  the  equations  by  the  quantity  K^/ViTa  ,  giving 
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The  out-of-plane  extension  of  a  pre-existing 
fracture  in  an  explosion-generated  stress  field.  The 
wing-cracks  have  grown  fo  a  length  I. 


a 


with 


S  j =  Cj  S3  "  Sq 


(3) 


S j  =  a,  fnZ/Kfc 
S3  =  a3  tfjTa/KIc 


(4) 


Crack. -GiQwtlL.and-InigriiciiQn 

Once  initiated,  the  wing  cracks  (as  we  shall  call  the  crack-like 
extensions  of  the  original  flaw)  grow  longer.  During  growth,  the 
stress  intensity  Ki  at  the  tip  of  each  wing  crack  is  equal  to,  or 
exceeds,  the  fracture  toughness  Kfc  of  the  solid.  The  condition  for 
crack  advance  is  simply 

Ki  *  Klc 

The  difference  between  tension  and  compression  is  that  growth  in 
compression  is  stable:  each  increment  of  crack  advance  requires  an 
increment  of  load,  at  least  until  the  cracks  start  to  interact  strongly. 

For  the  3-D  case  in  which  the  flaws  and  wing  cracks  are 
entirely  contained,  within  the  material,  Ashby  and  Sammis  (1989) 
calculated  Kj,  at  the  tip  of  the  wing  cracks  to  arrive  at  the  following 
expression  which  relates  the  length  of  the  wir.g  cracks  to  the  applied 
stresses 
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TABLE  1.  MATERIAL  PROPERTIES  OF  THE  ROCKS  AND  FITTED  MODEL  PARAMETERS  (Bracketted  data  are  estimates}. 
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(5) 


where  Sj  and  S3  arc  defined  by  equation  (4).  The  values  of  the 
constants  arc: 


c  (  i+h2)>/2  +  h 

( l+p2),'2-p 

C2  =  na  (( 1 4-  ja2) U2-  p) 

C3  =  2 

c4  = 2  (( 1 +^2) m-  J 


(6) 


p  =  0.1 


The  preexisting  damage,  D0,  and  the  current  damage,  D,  are  defined 
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as 


D0  =  i7c(aa)3Nv  (7) 

D»j;t(l  +  aa)3Nv 

where  I  is  the  length  of  the  growing  wing-cracks  and  Nv  is  the 
number  of  initial  cracks  of  length  a  per  unit  volume. 

Figures  3  and  4  illustrate  some  features  of  the  resu-s.  The 
axial  stress  at  first  rises  as  damage  grows  (Figure  3),  passing  through 
a  peak  which  shifts  to  higher  values  of  damage  as  the  confining 
pressure  increases.  The  damage  surface,  shown  in  Figure  4,  is  almost 
l  cone  meaning  that,  to  a  first  approximation,  the  failure  envelope  is 
described  by 


°>  -  Co>  -  °c  (8) 

where  C  is  a  constant  and  <Jc  is  the  unconfined  compressive  strength. 
The  value  of  the  model  is  that  it  gives  a  physical  interpretation  to  C 
and  0c  »  and  relates  them  to  the  initial  uamar ,,  me  coefficient  of 
friction,  and  the  crack  size. 

Analysis-  of  -Pata 

The  strength  of  many  different  rocks  have  been  determined 
under  triaxial  loading  conditions  (  loll  >  I02I  =  I03I,  where  all  three 
principal  stresses  are  compressive)  .  For  a  few,  the  initiation  of 
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Figure  3.  The  dependence  of  axial  stress  on  damage  a 
predicted  by  equations  (5)  -  (7).  The  peak  stress  is 
marked.  We  take  this  as  the  failure  stress. 
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Figure  4.  The  peak  value  of  axial  stress  plotted  against  the 
confining  stress  to  show  the  failure  surface.  In  this 
formulation,  damage  is  a  state  variable. 
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microcracking  has  also  been  determined.  Ashby  and  Sammis  (1989) 
applied  the  damage  mechanics  model  developed  above  to  nine 
different  rock  types  for  which  the  most  complete  data  sets  exist  : 
granite,  aplite,  duniic,  eclogitc,  gabbro,  sandstone,  limestone,  marble 
and  rock  salt.  These  rocks  represent  a  wide  range  of  composition 
and  initial  damage.  Granite,  aplite,  dunite  and  eclogitc  arc  low 
porosity,  crystalline,  igneous  rocks  in  which  the  initial  damage  is 
mostly  in  the  form  of  low  aspect  cracks.  Limestone,  and  sandstone 
arc  porous  sedimentary  rocks  in  which  the  initial  damage  is  mostly 
in  the  form  of  high  aspect  pores.  Marble  is  a  mctamorphic  rock  with 
initial  damage  of  a  form  intermediate  between  the  previous  two 
extremes.  These  rocks  also  span  a  wide  range  of  yield  strengths. 

The  igneous  rocks  have  yield  strengths  in  excess  of  2  GPa  while  the 
calcareous  rocks  (limestone  and  marble)  yield  at  stresses  below  1 
GPa.  Rock  salt  is  at  the  low  extreme  with  a  yield  strength  below  100 
MPa. 

The  damage  mechanics  model  formulated  above  has  three 
constants:  p  t  a,  and  (1.  They  arc  not  strictly  adjustable  because  the 
theory  makes  predictions  for  their  values.  However  because  of 
approximations  in  the  derivation,  and  the  uncertainties  in  the  aspect 
ratio  of  the  starter  flaws,  we  have  treated  p  as  adjustable,  choosing 
the  value  0.45  to  give  the  best  fit  to  the  data. 

For  those  materials  where  crack  initiation  data  are  available, 
the  crack  length  2a  and  the  coefficient  of  friction  p  are  determined 
from  the  initiation  surface  using  eqn.  (2).  Klc  is  also  required. 
Although  Kic  may  be  estimated  for  most  rocks  (Atkinson  and 
Meredith,  1987),  the  starter  flaw  size  is  not  known  in  most  cases 


13 


and  must  bo  treated  as  an  adjustable  parameter.  The  derivation  of 
the  fundamental  equations  gives  C3  =  2,  and  this  gives  a  good 
description  of  the  materials  we  have  examined. 

As  an  example,  the  triaxial  data  for  damage  initiation  and 
failure  in  Westerly  granite  are  presented  on  plots  of  oi  vs  03  in 
Figure  5.  The  theoretical  fracture  initiation  surface  (eqn.  2),  surfaces 
of  constant  damage,  and  the  failure  surface  (calculated  from  the 
maximum  of  eqn.  (5)  for  each  value  of  03  )  arc  plotted  on  these 
graphs  for  comparison  with  the  data. 

As  the  confining  pressure  is  increased,  brittle  fracture  is  made 
increasingly  difficult.  A  critical  pressure  may  be  reached  at  which 
true  plasticity  replaces  crack  extension.  This  transition  can  be 
illustrated  by  plotting  a  yield  (or  creep)  surface,  defined  by: 

2  ,  f  2  2  2"] 

-  2  L(a|-a2>  +(a2-a3)  +(<Vcr,)J  (9) 

The  yield  surface  is  plotted  as  a  pair  of  heavy  broken  lines  on  the 
figure  The  yield  strength,  ay,  can  be  derived  from  hardness,  H,  data 
since  ay  =  H/3.  Figure  5a  shows  theoretical  surfaces  for  initiation, 
constant  damage,  oo.d  failure  at  low  values  of  the  confining  stress. 

The  fracture  initiation  data  from  Brace  et  al.  (1966)  were  determined 
from  the  onset  of  non-linear  behavior  of  the  volume  strain.  Only 
data  taken  at  the  highest  loading  rate  are  plotted  here  in  order  to 
minimize  effects  of  subcritical  crack  growth  (which  we  do  not  model). 
The  three  initiation  points  from  Holcomb  and  Costin  (1986)  were 
determined  from  the  onset  of  acoustic  emission  (AE)  in  a  previously 
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AXIAL  STRESS,  o',  (MPa) 


NORMALIZED  AXIAL  STRESS,  or,  /Kxl0"3 


Figure  5a.  Comparison  between  experimental  and  theoretical 
failure  surfaces  for  granite  at  low  and  intermediate 
confining  stress.  Data  and  theory  for  microfracture 
initiation  and  surfaces  of  constant  damage  are  also 
compared. 
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RADIAL  STRESS,  cr3  (MPa) 


NORMALIZED  RADIAL  STRESS  o$/KxKf 


AXIAL  STRESS  o’,  (MPa) 


NORMALIZED  AXIAL  STRESS  ajK  xIO"3 


Figure  5b.  Comparison  between  experimental  and  theoretical 
failure  surfaces  for  granite  exteded  to  the  largest 
measured  confining  stress.  The  yield  surface  (equation  9) 
s  also  plotted  as  the  heavy  broken  *ines. 
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RADIAL  STRESS  aj  (MPa) 


unstressed  sample.  Also  shown  is  a  surface  of  constant  damage 
mapped  by  Holcomb  and  Costin  (1986)  using  the  AE  Kaiser  effect  as  a 
probe.  The  triaxifd  failure  data  arc  from  Brace  et  al..  (1966)  and 
Mogi  (1966). 

Figure  5b  shows  the  failure  surface  extended  out  to  large 
values  of  the  confining  stress  where  it  intersects  the  yield  surface. 
Data  at  low  and  intermediate  confining  pressures  arc  from  the  same 
sources  as  in  Figure  5a.  Those  at  high  confining  pressures  arc  from 
Schock  and  Heard  (1974)  for  Westerly  granite  and  from  Shimada 
(1981)  for  Man-nari  granite  (grain  size  1-3  mm,  apparent  porosity: 
0.7%).  It  is  evident  that  the  failure  surface  has  considerable 
curvature  and  deviates  from  our  theoretical  model  at  high  confining 
pressures.  Although  Janach  and  Guex  (1980)  have  modeled  this 
curvature  in  terms  of  the  formation  of  shear  bubbles  at  the  grain 
boundaries,  Figure  5b  supports  the  possibility  that  the  curvature  is 
due  to  a  gradual  transition  to  ductile  behavior.  Analogous  curvature 
is  evident  in  figures  for  limestone,  marble,  and  NaCl  (Ashby  and 
Sammis,  1989)  which  are  known  to  exhibit  ductile  behavior  at 
moderate  confining  pressure,  although  the  curvature  in  these  rocks 
occurs  over  a  more  limited  pressure  range.  The  broader  transition  in 
granite  may  reflect  its  multimincralic  composition  for  which  the 
individual  minerals  have  different  brittle  -  ductile  transition 
pressures.  Note  that  the  two  granites  show  different  transitional 
behavior  at  intermediate  pressures,  but  that  both  approach  the  same 
ductile  limit.  Schock  and  Heard's  observation  that  the  stress-  strain 
curve  is  linear  to  failure  may  be  due  to  the  convergence  of  the 
initiation  and  failure  surfaces  at  very  high  confining  pressures  in 
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Figure  5b.  The  sudden  release  of  energy  and  shear  localization  at 
failure  do  not  preclude  stress  concentration  by  ductile  processes  in 
the  weaker  minerals. 

AP-PLICATiQN-TO  LABORATORY  SIMULATORS  JM.GRAN1TB 

The  damage  mechanics  approach  may  be  used  to  interpret  the 
results  of  a  recent  set  of  laboratory  simulations  in  granite  which 
were  designed  to  investigate  the  pulse-broadening  observed  in  the 
field  (Nagy  and  Florence,  1986).  The  results  of  these  experiments 
arc  stmmarized  in  Table  2. 


TABLE  2 :  LABORATORY  SIMULATIONS 
Sample  Description  EuisC-Shapg 


As  received  (dry) .  Narrow 

(wet) .  Narrow 

Thermally  fractured  (dry) . . .  Narrow 

(wet) .  Narrow 

Gas  fractured  (dry) .  Narrow 

(wet) .  Slightly  wider 

(2-3x) 


The  gas  fractured  samples  were  the  most  extensively  fractured, 
and  probably  had  the  largest  fractures  since  they  were  barely 
cohesive.  Only  these  samples  which  contained  water  showed  any 
pulse  broadening.  All  other  cases  produced  a  narrow  pulse  in 
agreement  with  numerical  predictions  based  on  strength  parameters 
measured  in  the  laboratory. 


1  8 


Wc  interpret  the  lack  of  pulse  broadening  in  these  simulations 
to  be  the  consequence  of  a  scaling  effect  which  is  inherent  in  the 
physics  of  crack  growth.  The  stress  necessary  to  nucleate  and  extend 
out-of-plane  crack  growth  scales  as  the  inverse  square  root  of  the 
initial  flaw  size.  Since  the  maximum  initial  flaw  size  in  the 
simulations  is  limited  by  the  sample  dimensions,  the  experimental 
results  can  be  explained  if  these  flaws  arc  too  short  to  be  activated  at 
the  stress  levels  in  the  laboratory  explosion. 

Although  quantitative  verification  of  this  explanation  will 
require  a  full  numerical  simulation  (which  is  in  progress),  the 
following  approximate  calculation  illustrates  the  point.  Consider  the 
stress  pulse  at  a  distance  of  204  m  from  the  Piledriver  explosion 
shown  in  Figure  6.  This  pulse  was  calculated  by  Rimer  et  al.  (1987) 
using  laboratory  strength  data  for  granite.  At  the  stress  peak,  crrr  = 
334  MPa  and  aoo  =  135  MPa.  Equations  (5)  and  (6)  were  used  to 
calculate  the  crack  growth  under  this  str.ss  state  and  the  results  are 
shown  in  Figure  7a.  For  the  dry  case,  it  is  expected  that  p=0.6.  For 
the  saturated  case,  it  is  assumed  that  the  water  will  hold  the  cracks 
open  and  p  =  0  .  For  either  value  of  p,  the  hoop  stress  is  large  enough 
that  significant  crack  growth  does  not  occur. 

After  the  peak,  the  stresses  fall  to  values  of  orr  =  33  MPa  and 
009  =  2  MPa.  Crack  growth  under  these  conditions  is  shown  in  Figure 
7b.  Note  that  in  the  dry  case,  cracks  larger  than  about  one  meter 
show  significant  growth.  Hence  a  model  simulation  using  a  sample 
with  dimensions  approximately  equal  to  or  smaller  than  one  meter 
would  not  show  weakening.  For  the  wet  case,  cracks  larger  than 
about  one  centimeter  show  significant  growth.  Hence  a  meter  size 
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sample  might  be  expected  to  show  some  weakening  effects  if  it  were 
saturated.  This  is  in  agreement  with  the  experimental  observations 
summarized  in  Table  2. 

SEISMIC  PULSE  BROADENING  ASSOCIATED  WITH  FRACTURE  DAMAGE 
Equation  (5)  relates  the  current  damage  D  to  the  state  of  stress. 
It  can  be  used  to  generate  a  stress-strain  curve  if  the  elastic  modulus 
is  also  known  as  a  function  of  the  damage.  Saminis  and  Ashby 
(1986)  show  how  the  tangent  moduli  can  be  evaluated  from  the 
stress  i».>jnsity  factor  to  yield  the  following  expressions  for  the 
radial  modulus  Er  and  the  axial  modulus  E»  : 


E,-Eo(l-D) 


(10) 


E,«Eo 


-  jtLrl 

I 


(ii) 


where 

Stress-strain  curves  calculated  using  equations  (5)  and  (10)  - 
(12)  have  the  general  form  shown  in  Figure  8c.  As  the  rock  is 
loaded,  the  microcracks  extend  out-of-plane  as  discussed  above,  that 
is,  their  length  increases  from  Lo  to  Ll  Crack  growth  and 
interaction  produce  a  nonlinear  stress-strain  path.  When  the  rock  is 
unloaded,  growth  stops,  and  the  crack  length  remains  at  Lj  during 
the  unloading  process.  If  the  rock  as  a  whole  still  behaves  as  an 
elastic  body,  then  we  have  oi=  E(Li)e  for  the  unloading  process. 
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stresses  (kbar) 


Figure  6.  Stress  pulse  at  a  distance  of  204  m.  from  the 

Piledriver  explosion  calculated  by  Rimer  et  al.  (1987) 
using  laboratory  strength  parameters  for  granite. 
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Crack  Extension,  i/a  C:  Crack  Extension, 


Raw  Size,  a  (mm) 


Flaw  Size,  a  (mm) 

Figure  7.  Crack  extension  as  a  function  of  the  initial  crack 

size.  Case  (a)  is  at  the  peak  stress  in  Fig.  6  while  case  (b) 
is  at  a  post-peak  stress  state. 
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Figure  8.  Stress  strain  curves  and  pulse-shape  distortion  for 

(a)  the  linear  elastic  case,  (b)  the  non-linear  elastic  case, 
and  (c)  the  damage  mechanics  case.  Note  that  only  the 
damage  case  produces  pulse  broadening. 


Noia  that  the  unload  is  controlled  by  the  secant  modulus,  not  the 
tangent  modulus  which  controls  loading. 

In  order  to  investigate  the  effect  of  a  damage-based  rheology 
on  the  waveform,  we  assume  a  stress-strain  curve  which  has  the 
above  behavior  and  is  of  the  form: 


cr*Eo(l -e-*rt«)  for  loading  (13) 

<j >=F.q(i for  unloading 


The  basic  equations  which  govern  wave  propagation  in  a 
damaged  medium  ar  the  equations  of  motion: 


92u  do  P, 


(14) 


The  difference  equations  corresponding  to  eqns.  (10)  and  (11)  are: 


=  2uh  -  +  fa2  Ax  otn+j^2  -  ojj.j/J  +  P„ 


(15) 


°jwl/2  ®C  £^+i/2 
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This  formulation  is  a  general  one  which  includes  the  linear  clastic 
case,  the  non-linear  clastic  ease,  and  the  damage  case.  The  linear 
case  was  tested  against  the  known  analytic  solution  for  the  case 
where  i(x,  t  )  is  a  Gaussian  function. 

Waveforms  calculated  for  the  clastic,  non-lir.car  clastic,  and 
damage  eases  arc  compared  in  Figure  8.  For  the  clastic  ease,  the 
shape  of  the  response  pulse  is  the  same  as  that  of  the  input  force.  As 
the  input  force  increases  for  0  to  co,  the  response  increases  from 
point  1  to  point  2.  As  the  input  force  decreases  from  gq  to  0,  the 
response  is  from  point  2  to  3.  For  the  damage  case,  the  response 
pulse  can  be  divided  into  three  segments.  For  the  first  segment, 
from  point  1  to  2,  the  response  is  almost  the  same  as  the  linear 
elastic  case.  The  second  segment,  from  point  2  to  3,  deviates  from 
the  elastic  ease,  in  the  sense  that  it  propagates  more  slowly  due  to 
the  change  in  slope  of  the  stress-strain  curve.  The  last  segment,  the 
unloading  from  3  to  4  is  clastic  where  the  modulus  is  the  secant 
modulus  shown  in  Fig.  3.  The  result  is  that  the  response  pulse  for 
the  damage  case  is  wider  than  that  for  the  elastic  or  the  non-linear 
elastic  cases. 

PiS-CUSSlQN 

We  have  developed  a  damage  mechanics  for  compressive  stress 
which  has  been  verified  using  existing  data  for  a  wide  variety  of 
rocks  (Ashby  and  Sammis,  1989).  A  simple  one-dimensional 
application  of  this  model  provides  a  plausible  explanation  of  the 
pulse  broadening  observed  for  underground  explosions  in  granite. 

The  scaling  relation  between  the  size  of  the  preexisting  flaws  and  the 
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stress  required  to  increase  damage  offers  an  explanation  of  why  such 
broadening  is  neither  predicted  by  the  rheological  parameters 
measured  on  granite  specimens  in  the  laboratory  nor  observed  in 
laboratory  explosion  simulations. 
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ABSTRACT 

The  development  of  microcrack  damage  in  brittle  solids  in  compression 
is  analyzed,  using  a  simple  model.  The  model  is  developed  from  recent  detailed 
analysis  of  the  initialtion,  propagation  and  linkage  of  microfracturcs  from 
pre-existing  cracks,  voids,  or  other  inhomogcnctics.  It  describes  the  evolution 
of  damage  with  strain  and  from  it  the  stress-strain  response,  and  criteria  for 
failure,  can  be  established.  The  results  arc  used  to  construct  failure  surfaces 
in  stress  space  which  combine  information  about  brittle  failure  with  data 
describing  the  onset  of  plastic  yielding.  Such  failure  surfaces  arc  constructed 
for  a  number  of  rocks  and  are  compared  with  previously  published 
experimental  data. 
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1.  INTRODUCTION 


When  a  brittle  .solid  is  loaded  (o  failure,  it  docs  so  by  the  propagation  of 
cracks.  The  cracks  nucleate  and  propagate  from  inhomogcncitics.  by  which  we 
mean  holes,  inclusions,  microcracks,  surface  scratches  or  other  defects.  The 
difference  between  compressive  and  tensile  fracture  is  that  in  tension  a  single 
crack  grows  unstably  (once  started,  it  accelerates  across  the  sample  to  cause 
failure)  while  in  compression  a  population  of  small  cracks  extend  stably,  each 
growing  longer  as  the  stress  is  raised,  until  they  interact  in  some  cooperative  way 
to  give  final  failure  (Figure  1).  Because  of  this,  the  strength  of  a  brittle  solid  in 
compression  is  usually  greater,  by  a  factor  of  ten  or  more,  than  that  in  tension. 

Measurements  of  the  crushing  strength  of  stone,  brick  and  of  cement  must 
have  been  of  interest  to  civil  engineers  since  pre-Roman  times.  Systematic 
measurements  of  compressive  strength  really  began  about  the  middle  of  the  last 
century  (for  its  history,  see  Jaeger  and  Cook,  1976)  but  without  much  attempt  to 
understand  what  determined  it,  or  why  brittle  materials  had  useful  strength  in 
compression  but  none  to  speak  of  in  tension.  Elucidation  of  the  mechanics  of 
brittle  tensile  fracture  has  its  roots  in  the  work  of  Griffith  (1924),  Irwin(1958) 
and  others  that  followed  (sec  Knott,  1973,  for  a  review),  which  has  led  to  the 
development  of  fracture  mechanics  as  a  branch  of  engineering  design.  The 
understanding  of  compressive  brittle  fracture  is  more  recent,  and  still 
incomplete.  A  recent  series  of  papers  and  reviews  (Griggs  and  Handin,  1960; 
Paterson,  1978;  Ilallbaucr  el  al„  1973;  Tapponnier  and  Brace,  1976;  Wawersik  and 
Fairhurst,  1970;  Wawersik  and  Brace.  1971;  Ncmat-Nasscr  and  Horii,  1982;  Newman. 
1978;  Ashby  and  Hallam,  1986;  Sammis  and  Ashby,  1986)  have  established  that  an 
isolated  crack  in  a  large  body  grows  stably  until  its  length  becomes  comparable 
with  the  dimension  of  the  body  itself;  and  that  when  many  cracks  arc  present  (as 
they  always  arc  in  natural  rocks,  in  brick,  in  concrete  and  most  ceramics)  the 
cracks  grow  stably  until  their  length  is  comparable  to  their  spacing  when  they 
interact,  an  instability  develops,  and  the  sample  fails. 
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The  problem  can  be  complicated  by  time-dependent  effects  (Anderson  and 
Grew,  1976;  Martin,  1972;  Waza  ct  al„  1980;  Sano  ct  al„  1981,  Costin  and  Holcomb, 
1981;  Atkinson  and  Meredith,  1987b),  which  have  at  least  two  origins.  On  the  one 
hand  crack  growth  can  be  limited  by  a  chemical  reaction,  often  with  water.  On 
the  other,  cracking  in  compression  is  associated  with  dilation;  if  the  body  is 
saturated  with  a  fluid,  then  its  flow  into  the  dilating  region  can  introduce  a  time- 
dependent  aspect  to  fracture.  In  both  eases,  a  static  load  which  docs  not 
immediately  cause  failure  may  still  do  so  if  left  in  place  for  a  sufficient  length  of 
lime. 

The  understanding  of  compressive  brittle  fracture  is  still  incomplete,  but 
the  mechanisms  involved  arc  much  clearer  than  a  decade  ago.  It  seems  an 
appropriate  time  to  try  to  abstract  from  the  new  observations  and  modelling  a 
simplified  description  of  compression-cracking,  basing  it  as  far  as  possible  on  the 
physical  understanding.  The  goal  is  to  develop  a  damage  mechanics  of  brittle 
solids,  from  which  the  stress-strain  response  and  an  operational  definition  of 
failure  can  be  derived  for  a  material  with  a  given  set  of  elastic  properties  and 
given  defect  population,  under  a  given  state  of  stress.  Two  attempts  to  achieve 
this  can  be  found  in  the  open  literature;  that  of  Costin  (1983,  1985).  and  that  of 
Sammis  and  Ashby  (1986).  Central  to  the  problem  is  the  relationship  between 
stress  and  crack  extension.  Costin  (1983,  1985)  postulates  a  relationship  of 
reasonable  form,  and  develops  from  it  expressions  for  the  failure  surface  which 
(with  some  adjustable  parameters)  give  a  good  description  of  the  experimental 
data  available  at  that  time,  but  the  model  is  not  based  on  a  physical  model  for 
crack  growth.  Sammis  and  Ashby  (1986)  and  Ashby  and  Hallam  (1986)  use 
methods  of  fracture  mechanics  to  develop  a  physical  model  for  crack  extension, 
which  they  use  to  plot  stress-strain  curves  for  brittle  solids  from  which  failure 
surfaces  can  be  constructed,  but  the  complexity  of  their  model  makes  the  process 
cumbersome.  In  the  present  paper  we  attempt  to  develop  a  simpler,  model-based 
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mechanics  of  brittle  compressive  fracture,  drawing  heavily  on  the  previous 
pieces  of  work. 

2.  CRACK  INITIATION. IN  CQMERESSIQK 

Most  brittle  solids  contain  inhomogcncitics:  small  holes  or  cracks,  particles 
which  arc  poorly  bonded,  or  phases  which  have  different  moduli  or  strengths 
from  those  of  the  matrix.  Any  one  of  these  can  act  as  nuclei  for  new  cracks  when 
the  solid  is  loaded. 

The  range  of  possible  nuclei  is  wide,  but  the  spectrum  of  their 
characteristics  is  probably  bracketed  by  two  extremes:  the  spherical  hole  and  the 
sharp  inclined  crack  (Figure  2).  Both  have  been  studied  experimentally  and  both 
have  been  modeled,  the  first  by  Sammis  and  Ashby  (1986)  and  the  second  by 
Nemat-Nasscr  and  Horii  (1982)  and  Ashby  and  Hallam  (1986).  In  both  eases,  the 
criterion  for  crack  initiation,  under  axisymmetric  loading  has  the  form 

CTI  =  C1  <*3  "  a0 

where  ci  and  ero  an  material  properties,  oj  is  the  axial  stress,  and  02*03  the  radial 
stress  (both  positive  when  tensile,  negative  when  compressive). 

In  the  later  development  of  this  paper  we  consider  the  growth  of  crack- 
damage  from  initial,  inclined  cracks  as  in  Figure  2a.  For  this  case  (Nemat-Nasser 
and  Horii,  1982;  Ashby  and  Hallam,  1986)  cracks  initiate  when 


where  it  is  the  coefficient  of  friction  acting  across  the  crack  faces,  Kic  is  the 
fracture  toughness  of  the  material  through  which  the  new  crack  propagates,  and 
2a  is  the  length  of  the  original  inclined  crack.  Rocks,  typically,  show  a 
coefficient  of  friction  of  about  0.6,  in  which  ease  Cl  =3.1  and  cr0=3.1  Kk/ Vrc  a  . 

Crack  initiation  from  holes  (Sammis  and  Ashby,  1986)  gives  similar  values. 
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Crack  initiation  can  be  detected  in  several  ways:  by  the  start  of 
acoustic  emission,  by  the  first  non-linearity  of  the  stress-strain  curve,  by  the 
dilation  of  the  sample,  or  by  a  sudden  increase  in  internal  friction.  None  give 
very  accurate  data,  but  they  do  allow  a  test  of  cqn.  (I).  Figure  3  shows  data  for 
crack  initiation  in  Westerly  granite  obtained  by  the  first  three  techniques 
(Holcomb  and  Costin,  19S6;  Brace  ct  al.,  1966)  plotted  on  axes  of  oj  and  03  to  allow 
comparison  with  cqn.  (1).  The  linear  relationship  gives  a  good  description  of  the 
data  with  a  slope  between  2.7  and  3.3  (corresponding  to  p  =  0.55  to  0.64)  and  an 
intercept  of  70  -  79  MPa  (corresponding  to  a  crack  length  2a  close  to  1  mm  when 
Kjc  =  1  MPa  mU2  ). 

The  theory  gives  an  adequate  description  of  the  data.  It  is  used  to  describe 
the  initiation  of  damage  in  the  diagrams  shown  later.  In  each  case,  experimental 
data  arc  fitted  to  cqn  '1)  to  give  p  and  a  (using  published  data  for  Kjc).  Results  of 
this  analysis  arc  summarised  in  Table  2.  In  the  computations,  it  is  convenient  to 
normalize  the  equations  by  the  quantity  K^/ViTa  ,  giving 


with 


S  j =  Cj  S3  ■  Sq 


Si  =  a,  Vna/Kjj 
S3  =  ct3  V7ta/Kjc 


(2) 

(3) 


S 


o  = 


3.  CRACK  GROWTH  AND  INTER  ACTION 
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Once  initialed,  the  wing  cracks  (as  we  shall  call  (he  crack-like  extensions 
of  ihc  original  flaw)  grow  longer.  During  growth,  the  stress  intensity  Kt  at  the 
tip  of  each  wing  crack  is  equal  to,  or  exceeds,  the  fracture  toughness  Kic  of  the 
solid.  The  condition  for  crack  advance  is  simply 

Kli  Kic 

The  difference  between  tension  and  compression,  as  already  mentioned,  is  that 
growth  in  compression  is  stable:  each  increment  of  crack  advance  requires  an 
increment  of  load,  at  least  until  the  cracks  start  to  interact  strongly.  We  will 
assume  that  a  steadily  increasing  load  drives  the  cracks  at  a  steady  rate,  though  in 
reality  the  inhomogcncily  of  natural  materials  may  cause  them  to  extend  in  little 
jumps.  The  problem,  then,  is  to  calculate  Kj  at  the  tip  of  the  wing  cracks. 
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3.1  Crack-Growth  -from  Starter  Flaws,  of  a  Single  Si/c  :  The  2-Dtmcnsional  Ca<e 
Figure  4  shows  an  array  of  through-cracks,  growing  in  a  linear-elastic 
medium  under  a  iriaxial  stress  field  a  ]  .  03  ,  positive  when  tensile,  negative  when 
compressive.  Consider  fitst  the  growth  of  a  single,  isolated  crack  from  an  initial 
inclined  flaw;  interaction  comes  later.  The  upper  inset  of  Figure  4  isolates  one 
crack:  it  is  made  up  ^f  an  initial  crack  of  length  2a  lying  at  an  angle  y  to  the  X( 
direction  with  two  wings,  each  of  length  I  which  (we  will  assume)  lie  parallel  to 
X  1 .  The  stress  intensity  at  '.he  tips  of  the  wing.;  is  obtained  approximately,  but 
adequately,  in  the  following  way,  based  on  the  work  of  Ncmat-Nasscr  and  Horii 
(1982),  Ashby  and  Hallam  (1986),  Horii  arj  Ncmat-Nasscr  (1985,  1986a  and 
Kemeny  and  Cook  (1986). 

'I  he  remote  field  01,03  creates  a  shear  stress  x  and  a  normal  stress  a  across 
faces  of  the  initial  crack.  The  crack  slides  (resisted  by  the  coefficient  of  friction 
p).  wedging  open  the  mouth  of  each  wing  crack  by  8  (Figure  4).  The  wedging  can 
be  thought  of  as  caused  by  forces,  F3,  parallel  to  X3,  acting  at  the  midpoint  of  the 
crack.  The  stresses  t  and  a  arc  given  by 

°3“cr, 

x  =  __sin2v  (4a) 


On  +  0,  0j—  0, 

0= - + - cos  2\p 

2  2 

F3  is  simply  the  component  of  the  sliding  force  acting  parallel  to  X3: 


(4b) 


0  r 


F3=  (T  +  p0)  2a  sin\j/ 
F3  =  -  (Aj  0!  -  A303)  a 


(5a) 


(5b) 


where  Ai  and  A3  arc  constants  which  depend  on  y,  to  be  determined  in  a  moment. 
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The  force  F3  ccting  ai  (he  mid  point  of  a  crack  of  length  21  creates  a  stress 


intensity  tending  to  open  the  crack  (Tada  et  ai.,  1985,  page  5.1)  of 


ViTT 


This  result  gives  a  good  estimate  of  the  stress  intensity  at  the  tip  of  a  wing  crack 
when  I  is  large,  but  it  breaks  down  (becoming  infinite)  when  I  is  vanishingly 
small.  The  stress  intensity  at  the  tip  of  the  initial  inclined  crack  is  not  infinite, 
but  can  be  calculated  exactly  as  explained  in  the  last  section.  We  overcome  this 
problem  by  introducing  an  "effective"  crack  length  (  I  +  pa)  giving 


Vic  (1  + pi) 


(6) 


We  then  choose  p  so  that  (Kj)i  becomes  equal  to  that  for  the  inclined  crack  when  1 
is  zero. 

Before  doing  this,  we  note  that  there  is  another  contribution  to  K[  at  the  tip 
of  the  wing  crack.  The  remote  confining  stress  03  acts  not  just  on  the  angled 
crack  but  on  the  wing  cracks  of  length  1.  In  so  doing,  it  produces  an  additional 
contribution  to  the  stress  intensity,  tending  to  close  the  crack  when  03  is 
compressive  (Tada  et  a!.,  1985,  page  5.1): 

(K|)3=o3V7T  (7) 


Summing  the  two  contributions,  with  Fg  given  by  cqn  (5b),  gives: 


K,  =.  —  ?■  — +  a3Vjt7 
Vrc(l  +  Pa) 


(8a) 
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(8b) 


A, a,  Vita  , —  f  a3  -A 

=  -  - ; . - . +  a^-YTta  - ~=r+  YL  | 

irvL+p  ^xvL-l-p  J 

where  Ls  I  /a.  The  eracks  extend  uniil  Kj  becomes  equal  10  Kjc. 

The  constants  arc  found  by  ensuring  that  this  equation  reduces  to  (he  exact 
result  for  crack  initiation  (L  =  0)  and  matches  the  known  results  for  very  long 
cracks  (L  »  1).  given  by  Nemat-Nasscr  and  Horii  (1982)  and  Ashby  and  Hallam 
(19S6),  cqns.  (3)  and  (6).  This  gives 

A  *^f(\  JV'2  ) 

A,=-=?-  VU+hJ  -w  (9, 


A3  =  Aj 


P  =  0.1 


Equation  (8b)  with  the  values  of  Aj,  A3,  and  P  (and  with  Kj  =  Kjc  )  is  plotted  in 
Figure  5.  It  shows  a,  VnTa  /Kk  plotted  against  L  with  the  earlier  numerical 

results  of  Nemat-Nasscr  and  Horii  (1982).  Equation  (8)  is  obviously  a  good 
approximation  to  the  earlier  calculations. 

Now  the  interaction.  The  main  part  of  Figure  4  shows  an  array  of  Na 
cracks  per  unit  area,  all  of  which  have  extended  to  a  length  2(  I  +  a  a.)  The  center- 
(o-ccntcr  spacing  of  the  cracks  is 


S  = 


JrA 


(10) 
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so  that  an  uncrackcd  ligament  of  average  length  S  *  2  (  1  +  oa)  remains  between 
the  cracks  in  the  X]  direction.  (Here  a  is  simply  a  geometric  constant,  and  must 
be  distinguished  from  P;  for  cracks  at  45°  to  Xi,  a»  1  /  V2* ) .  An  opening  force  F3 
acts  at  the  midpoint  of  each  crack.  Equilibrium  requires  that  this  opening  force 
be  balanced  by  a  mean  internal  stress  03  *  in  the  matrix,  as  shown  in  the  right 
hand  side  of  Figure  4.  The  average  internal  stress  is  given  by 


1 

°3 


rs 


S-2(l  +  att) 


CH) 


This  acts  on  the  wing  cracks,  so  that  cqn  (7)  now  becomes 


(12) 


We  now  define  the  initial  damage  Dq  and  the  current  damage  D  by: 


D0»x(aa)  Na 

(13a) 

D  =  jc(l+aa)2NA 

(13b) 

giving 


1/2 

t  -(AjCJ,-  A3o3)(D0/jO 

(14) 

a3“  (  l/2\ 

a  1 1  -2(D/7t)  ) 

Equations  (8a)  and  (8b)  now  become 
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Here  the  first  term  in  ihc  curly  brackets  describes  the  wedging  plus  the  crack- 
crack  interaction;  and  the  second  lean  describes  the  closing  effect  of  the  lateral 
confining  stress.  The  cracks  propagate  until  Ki  falls  to  Kic.  Using  this, 
rearranging  and  aggregating  the  constants  (with  V2/jt  «  1  )  gives  for 

proportional  loading  (with  X  =  03/01  held  constant)  : 


06) 


and  for  loading  at  constant  03 
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(17) 


where  S]  and  S3  arc  defined  by  equation  (3).  The  values  of  the  constants  arc 


(18) 


Figure  6  shows  how  the  axial  stress  cri  varies  with  damage  D  for  various  confining 
pressures.  The  left  hand  figure  shows  proportional  loading;  the  right  hand 
figure  ,  loading  at  constant  03.  The  peak  stress,  (ai)max  .  rises  and  moves  to  the 
right  as  X  or  <13  is  increased.  The  shapes  of  the  curves  at  constant  X  differ  from 


40 


those  at  constant  03,  as  expected,  but  the  peaks  arc  at  the  same  stress.  Figure  7 
shows  (aj)max  plotted  against  03  for  both  conditions:  the  points  lie  on  the  same 
line. 

3.2  Crack  Growth  from.  Starter  Flaws  of  a  Single  Size.  The  3-dimensional  Case 

It  is  usually  the  ease  that  flaws  arc  completely  contained  within  the 
material.  The  merit  of  the  2-dimcnsional  calculations  developed  in  section  3.2  is 
that  it  points  to  a  way  of  tackling  this  more  difficult  3-dimensional  problem.  We 
require  the  stress  intensity  at  the  periphery  of  a  contained  crack  emanating  from 
a  starter  flaw  (which  we  take  to  be  an  inclined,  penny-shaped  crack)  which  we 
will  equate,  as  before,  to  Kic  .  This  we  do  by  calculating  the  wedging  force  F3  as 
in  the  2-dimensional  ease.  The  wedging  force  creates  an  average  internal  stress 
03!  ,  The  stress  intensity  at  the  tip  of  a  given  wing  crack  is  calculated  from  the 
wedging  force  and  the  total  lateral  stress  (03  +  03*).  as  shown  in  Figure  8.  The 
significant  difference  in  the  results  of  this  3-D  and  the  earlier  2-D  calculations  is 
that  the  dependence  of  Si  and  S3  on  damage  D  involves  different  powers. 

The  wedging  force  F3,  as  before,  is  calculated  from  the  shear  and  normal 
stresses  (eqn  4)  acting  on  the  initial  crack  plane,  times  the  crack  area,  resolved 
into  the  X3  direction: 

F3=  (t+n<y)7ta2  siny 


=  -  (A,cr  A3cr3)a2  (19) 

Unlike  the  2-D  ease,  there  arc  no  exact  analytical  solutions  for  limiting  eases 
which  allow  Ai  and  A3  to  be  determined,  so  we  make  the  assumption  that  they 
have  the  same  values  as  before  (equation  9)  but  make  provision  to  adjust  them 
later  to  match  experimental  data. 

The  stress  intensity  Ki  has  the  same  three  contributions  as  before.  The 
wedging  force  F3  induces  a  stress  intensity  (Tada  ct  al„  1985,  page  24.2) 
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(20) 


(*0l 


J*3 


(tc(1  +  P»)) 


3/2 


where  j)  is  introduced  for  the  same  reason  as  before:  to  give  a  limiting  value  of 
(Ki> |  when  I  3  0.  The  contribution  due  to  03  and  03*  arc  (Tada  et  at.,  1983,  page 
24.2). 


(21) 


The  internal  stress  for  the  3*D  ease  is 


F, 


A-rt(l  +  cxa) 


(22) 


where  7t(J  +  aaJ  is  the  total  crack  area  projected  normal  to  X3  and  A  is  the  area 
per  crack 


A  = 


(23) 


where  Nv  is  the  number  of  cracks  per  unit  volume. 

Damage  is  defined  in  a  way  which  parallels  that  in  2-dimcnsions: 


4  3 

D0  =  j7t(aa)  Nv 


(24a) 


4  3 

D  =  j7t(!  +  aa)  Nv 


(24b) 


giving 
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(25) 


-  (AjOj- A3  (J3)Dq2/3 


The  stress  intensity  at  the  tip  of  the  wing  crack  is 


Gt(l  +  Pa)) 


3/2 


(26) 


2  3/2  f 

k  a  [(D/D0) 


1/3 


•  1  + 


1  +  2((D/D0),/3- 


2X 

A, 


As  before  the  cracks  propagate  until  Ki  falls  to  Kic.  Rearranging  and 
aggregating  the  constants  the  gives,  for  proportional  loading: 


and  for  loading  at  constant  a 3 
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(28) 


where  S]  and  S3  are  defined  by  equations  (3).  The  values  of  the  constants  arc, 


The  equations  and  constants  have  a  form  very  like  those  of  the  2- 
dimcnsional  model.  Two  significant  differences  should  be  noted.  First,  the  extra 
dimension  causes  the  powers  of  D  which  appear  in  the  equation  to  differ  (not 
surprisingly)  from  those  of  the  2-D  model.  Second,  the  constants  Cl  to  C4  arc  not 
known  with  the  same  precision  as  those  of  the  2-D  model  because  accurately- 
known  limiting  eases  are  not  available  to  calibrate  them.  We  shall  assume 
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(reasonably)  that  (he  dependence  on  the  coefficient,  4,  is  properly  included,  but 
that  the  constant  p  may  require  further  adjustment  to  give  a  good  match  with 
experiment. 

Figures  9  and  10  illustrate  some  features  of  the  results.  The  axial  stress  at 
first  rises  as  damage  grows  (Figure  9),  passing  through  a  peak  which  shifts  to 
higher  values  of  damage  as  the  confining  pressure  increases.  The  damage 
surface,  shown  in  Figure  10,  is  almost  a  cone  meaning  that,  to  a  first 
approximation,  the  failure  envelope  is  described  by 

a,  =  Ca3  -  ac 

where  C  is  a  constant  and  ac  is  the  unconfincd  compressive  strength.  The  value 
of  the  model  is  that  it  gives  a  physical  interpretation  to  C  and  oc  •  and  relates  them 
to  the  initial  damage,  the  coefficient  of  friction,  and  the  crack  size. 

4.  ANALYSIS  QF  DATA 

The  strength  of  many  different  rocks  have  been  determined  under  triaxial 
loading  conditions  (  lajl  >  la2l  =  I03I,  where  all  three  principal  stresses  arc 
compressive)  .  For  a  few,  the  initiation  of  microcracking  has  also  been 
determined.  We  now  apply  the  damage  mechanics  model  developed  above  to  nine 
different  rock  types  for  which  the  most  complete  data  sets  exist  :  granite,  aplitc, 
dunite,  cclogite,  gabbro,  sandstone,  limestone,  marble  and  rock  salt.  These  rocks 
represent  a  wide  range  of  composition  and  initial  damage.  Granite,  aplitc,  dunite 
and  cclogite  are  low  porosity,  crystalline,  igneous  rocks  in  which  the  initial 
damage  is  mostly  in  the  form  of  low  aspect  cracks.  Limestone  and  sandstone  arc 
porous  sedimentary  rocks  in  which  the  initial  damage  is  mostly  in  the  form  of 
high  aspect  pores.  Marble  is  a  metamorphic  rocks  with  initial  damage  of  a  form 
intermediate  between  the  previous  two  extremes.  These  rocks  also  span  a  wide 
range  of  yield  strength.  The  igneous  rocks  have  yield  strengths  in  excess  of  2  GPa 
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while  (he  calcareous  rocks  (limestone  and  marble)  yield  ni  stresses  below  1  GPa 
Rock  salt  is  at  the  low  extreme  with  a  yield  strength  below  100  MPa. 

The  damage  mechanics  model  formulated  above  has  three  constants:  p  ,  a, 
and  p.  They  arc  not  strictly  adjustable  because  the  theory  makes  predictions  for 
their  values.  However  because  of  approximations  in  the  derivation,  and  the 
uncertainties  in  the  aspect  ratio  of  the  starter  flaws,  we  have  treated  p  as 
adjustable,  choosing  the  value  0.45  to  give  the  best  fit  to  the  data. 

For  those  materials  where  crack  initiation  data  arc  available,  the  crack 
length  2a  and  the  coefficient  of  friction  p  are  determined  from  the  initiation 
surface  (sec  Table  2) 


ai  =  Cl  a3  +  o0 

where  Cl  and  o0  are  given  by  equation  (I).  Kjc  is  also  required.  Although  Kjc 
may  be  estimated  for  most  rocks  (Atkinson  and  Meredith,  1987a),  the  starter  flaw 
size  is  not  known  in  most  cases  and  must  be  treated  as  an  adjustable  parameter. 

The  derivation  of  the  fundamental  equations  (27)  and  (28)  from  equation  (26) 
gives  C3  =  2,  and  this  gives  a  good  description  of  the  materials  we  have  examined. 

Each  material  will  now  be  discussed  in  turn.  The  triaxial  data  for  damage 
initiation  and  failure  arc  presented  on  plots  of  01  vs  03  .  The  theoretical  fracture 
initiation  surface  (cqn.  1),  surfaces  of  constant  damage,  and  the  failure  surface 
(calculated  from  the  maximum  of  cqn.  (28)  for  each  value  of  03  )  arc  plotted  on 
these  graphs  for  comparison  with  the  data. 

As  the  confining  pressure  is  increased,  brittle  fracture  is  made 
increasingly  difficult.  A  critical  pressure  may  be  reached  at  which  true  plasticity 
replaces  crack  extension.  This  transition  can  be  illustrated  by  plotting  a  yield  (or 
creep)  surface,  defined  by: 
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(30) 


2  i  f  2  2  i] 

cfy  =  j|(ara2)  +  (a2-cr3)  +  (a3-a,)j 

The  yield  surface  is  ploucd  as  a  pair  of  heavy  broken  lines  on  each  figure  The 
yield  strength,  oy,  can  be  derived  front  hardness,  H,  data  since  cry  =  H/3.  The 
material  properties  and  constants  used  to  generate  the  theoretical  initiation  and 
failure  curves  arc  tabulated  for  each  solid. 

4.1  Granite 

Westerly  granite  is  a  fine-grained  (0.75  mm.),  low  porosity  (0.9%), 
isotropic,  two-mica  calc-alkaline  granite  which  has  become  a  standard  material  in 
rock  mechanics  testing  (see  Scholz,  1986,  for  a  b.'.cf  history).  Mincralogical 
modal  analyses  arc  given  by  Birch  (1960)  and  Wawersik  and  Brace  (1971). 

Figure  11a  shows  theoretical  surfaces  for  initiation,  constant  damage,  and 
failure  at  low  values  of  the  confining  stress.  The  fracture  initiation  data  from 
Brace  ct  al.  (1966)  were  determined  from  the  onset  of  non-linear  behavior  of  the 
volume  strain.  Only  data  taken  al  the  highest  loading  rate  arc  plotted  here  in 
order  to  minimize  effects  of  subcritical  crack  growth  (which  we  do  not  model). 

The  three  initiation  points  from  Holcomb  and  Costin  (1986)  were  determined  from 
the  onset  of  acoustic  emission  (AE)  in  a  previously  unstressed  sample.  Also  shown 
is  a  surface  of  constant  damage  mapped  by  Holcomb  and  Costin  (1986)  using  the  AE 
Kaiser  effect  as  a  probe.  The  triaxial  failure  data  arc  from  Brace  ct  al.  (1966)  and 
Mogi  (1966). 

Figure  lib  shows  the  failure  surface  extended  out  to  large  values  of  the 
confining  stress  where  it  intersects  the  yield  surface.  Data  at  low  and 
intermediate  confining  pressures  arc  from  the  same  sources  as  in  Figure  11a. 
Those  at  high  confining  pressures  arc  from  Schock  and  Heard  (1974)  for  Westerly 
granite  and  from  Shimada  (1981)  for  Man-nari  granite  (grain  size  1-3  mm, 
apparent  porosity:  0.7%).  It  is  evident  that  the  failure  surface  has  considerable 
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curvature  and  deviates  from  our  'heretical  model  at  high  confining  pressures. 
Although  Janach  and  Guex  (1980)  have  modeled  this  curvature  in  terms  of  the 
formation  of  shear  bubbles  at  the  grain  boundaries.  Figure  lib  supports  the 
possibility  that  the  curvature  is  due  to  a  gradual  transition  to  ductile  behavior. 
Analogous  curvature  is  evident  in  subsequent  figures  for  limestone,  raarblc,  and 
NaCI  which  arc  known  to  exhibit  ductile  behavior  at  moderate  confining 
pressure,  although  the  curvature  in  these  rocks  occurs  over  a  more  limited 
pressure  range.  The  broader  transition  in  granite  may  reflect  its  multimincralic 
composition  for  which  the  individual  minerals  have  different  brittle  -  ductile 
transition  pressures.  Note  that  the  two  granites  show  different  transitional 
behavior  at  intermediate  pressures,  but  that  both  approach  the  same  ductile  limit. 
Schock  and  Heard's  observation  that  the  stress-  strain  curve  is  linear  to  failure 
may  be  due  to  the  convergence  of  the  initiation  and  failure  surfaces  at  very  high 
confining  pressures  in  Figure  lib.  The  sudden  release  of  energy  and  shear 
localization  at  failure  do  not  preclude  stress  concentration  by  ductile  processes  in 
the  weaker  minerals. 

4.2  AjlLLIC 

Brace  et  al.  (1966)  studied  a  quartz-oligoclasc  aplitc  (63%  oligoclasc,  27% 
quartz,  and  10%  biotitc)  which  they  described  as  fine-grained  and  flinty, 
apparently  isotropic,  and  of  high  strength.  The  feldspar  is  highly  altered.  The 
grain  size  of  the  ground  mass  is  about  40  pm  and  of  the  phcnocrysts,  about  100  pm. 
The  small  grain  size  and  high  strength  arc  consistent  with  the  small  flaw  size 
required  by  the  model  (see  Figure  12).  The  relatively  high  initial  damage  is 
consistent  with  the  flinty  texture.  A  high  density  of  small  flaws  may  explain  why 
flinty  materials  can  be  rcl'ably  fashioned  into  tools  by  flaking  off  small  bits  in  a 
controllable  manner. 

4.3  Pun  its 
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Dunitc  is  an  almost  pure  olivine  rock.  Shimada  ct  al.  (1983)  measured  the 
compressive  strength  of  Horoman  dunitc  (grain  Size  Q.l-0.9  mm)  at  confining 
pressures  up  to  450  MPa  using  a  conventional  triaxial  testing  apparatus  and  to  3 
GPa  using  a  cubic  press.  Acoustic  emissions  showed  a  change  in  failure  mode  at 
confining  pressures  between  0.44  and  0.77  GPa.  Below  these  pressures,  AE  activity 
increased  rapidly  between  the  onset  of  dilatancy  and  failure  (the  typical  pattern 
for  brittle  failure).  Above  the  transition  pressure,  and  increase  in  AE  activity  was 
not  observed  to  precede  failure;  rather,  the  level  of  AE  remained  nearly  constant 
up  to  failure.  Shimada  ct  al.  (1983)  correlate  this  change  in  AE  behavior  with  the 
extreme  curvature  in  the  failure  envelope.  As  is  evident  in  Figure  13,  our  model 
suggest  that  this  change  in  behavior  is  associated  with  the  transition  to  plastic 
deformation. 

4.4  EdfiiilC 

Eclogitc  is  an  uliramafic  pyroxene-garnet  rock.  The  Akaishi  eclogitc 
measured  by  Shimada  ct  al.  (1983)  was  composed  of  0.1-0.3  mm  pyroxene  grains 
and  O.S-2.3  mm  garnet  grains.  It  had  a  density  of  3.642  gni/cm^  and  a  porosity  of 
0.4%.  Conventional  triaxial  tests  covered  a  range  of  confining  pressures  from  0  to 
450  MPa  while  tests  in  an  opposed  anvil  cubic  press  extended  the  confining 
pressure  to  3  GPa.  As  discussed  above  for  dunitc,  the  AE  patterns  indicate  a 
change  in  failure  mechanism  at  confining  pressures  between  1.02  and  1.99  GPa. 
The  data  and  theoretical  surfaces  arc  given  in  Figure  14. 

4.5  Oahhro 

In  addition  to  dunitc  and  eclogitc  discussed  above,  Shimada  ct  ai.  (1983)  also 
studied  Murotomisaki  gabbro,  a  hypcrsthcnc-bcaring-olivinc-augite  gabbro.  The 
grain  size  of  the  olivine  component  is  1-2  mm,  pyroxene  is  about  0.7  mm,  and  the 
plagioclasc  is  about  0.7-3  mm.  The  bulk  density  is  2.985  gm/cm^  and  the  reported 
porosity  is  0.4%.  The  pattern  of  AE  indicated  a  change  in  failure  mechanism  at 
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confining  pressures  between  0.5!  and  0.76  GPa.  The  data  and  theoretical  surfaces 
arc  given  in  Figure  15. 

4.6  SifidslQDfi 

The  sandstone  data  in  Figure  16  were  obtained  in  triaxial  compression  by 
Gowd  and  Rummct  (1980).  The  rock  is  described  as  a  medium  grain-sUed 
Bumsandstane  from  SW-Gcrmany  with  subangular  to  round  quart*  grains  bedded 
within  a  clayey  matrix..  Its  initial  porosity  was  15%  with  an  initial  permeability 
of  50  nticrodarcy.  The  damage  initiation  data  were  defined  by  the  onset  of 
dilatancy.  At  confining  pressures  above  about  30  MPa.  the  strcss*strain  curves 
are  noniincar  at  lower  values  of  the  axial  stress  than  the  observed  onset  of 
dilatancy.  This  probably  reflects  the  suppression  of  dilatancy  by  pore  collapse,  a 

phenomenon  which  the  authors  propose  to  explain  the  total  lack  of  observed 
dilatancy  at  the  highest  confining  pressures.  Such  effects  are  beyond  the  scope  of 
our  model. 

A  transition  from  brittle  failure  to  apparent  ductile  shear  deformation 
takes  place  at  a  pressure  of  about  100  MPa.  However,  the  observed  pressure 
dependence  of  the  flow  stress  for  confining  procures  in  excess  of  100  MPa  argues 
against  true  ductile  flow  and  for  a  cataclastic  mode  of  deformation  probably 
involving  pore  cc'dapsc.  Dilatancy  at  failure  is  a  constant  for  confining 
pressures  between  0  and  40  MPa.  From  40  to  100  MPa,  dilatancy  at  failure  decreases 
to  zero.  Above  100  MPa,  brittle  failure  does  not  occur. 

4.7  Limestone 

Solcnhofcn  limestone  is  a  fine  grained  (0.01  mm)  mechanically  isotropic 
limestone  from  Bavaria.  It  has  a  connected  porosity  of  5.3%  and  a  total  porosity  in 
the  range  6  -  9%  (Rutter,  1971’.).  The  strength  data  in  Figure  17  arc  from  Heard 
(1960)  and  include  both  triaxial  compression  and  tension.  The  fracture  initiation 
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points  were  picked  ns  the  onset  of  nonlinearity  in  his  published  stress-strain 
curves  and  are  only  approximate. 

*1.8  Maxhls 

The  only  marble  for  which  fracture  initiation  data  is  available  is  described 
by  Brace  ct  al.  (1966)  as  a  medium  grained  almost  pure  caleitc  marble  of  unknown 
origin.  They  report  it  to  be  apparently  isotropic,  very  ductile  even  at  low 
confining  pressures,  and  having  a  grain  size  of  about  0.2  mm. 

The  unusually  low  fraeturc  initiation  stress  (Figure  IS)  requires  either 
large  starter  flaws  u  6  mm)  or  a  low  fracture  toughness.  Since  Atkinson  and 
Meredith  (I987at  report  K(e  as  low  as  0.19  MPa  m^  for  caleitc,  we  have  fit  the 
initiation  data  using  this  value  which  then  implies  a  starter  flaw  size  of  0.4  mm, 
which  is  comparable  to  the  grain  size. 

In  a  microscopic  study  of  nuclcation  in  marble,  Olsson  and  Peng  (1976) 
found  that  microcracks  often  nucleate  where  slip  bands  intersect  grain 
boundaries.  Although  such  slip-bands  arc  physically  analogous  to  angle  cracks, 
there  may  be  a  significantly  larger  number  of  such  nuclei  since  every  favorably 
aligned  grain  is  a  potential  source  of  nuclei.  This  may  explain  the  large  values  of 
initial  damage  Do  required  to  fit  the  marble  data. 

The  data  set  which  we  fit  is  for  Carrara  marble  which  is  the  fine  grained** 
(about  0.1  mm)  isotropic  marble  used  by  Michaelangelo  for  the  Picta  and  other 
well-known  works.  Its  total  porosity  is  about  1.1%  (Edmond  and  Paterson,  1972). 
The  triaxial  data  in  Figure  18  are  from  Von  Karman  (1911)  and  Edmond  and 
Paterson  (1972).  Brittle  vjl  ductile  behavior  was  deduced  from  the  shape  of  the 
stress-strain  curve  and  from  the  volume  changes  associated  with  the 
deformation. 
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4.9  Rock-Salt 


Rock  sail  exhibits  a  room-temperature  briule-to-dueiile  transition  at  the 
lowest  confining  pressure  of  any  rock  in  this  study,  llunschc  (1981)  tested  three 
types  of  natural  salt  at  low  confining  pressure  under  both  the  common  triaxial 
loading  and  the  less  common  true  multiaxial  loading  at  strain  rates  of  about  10’®  s* 
K  llandin  (19531  collected  conventional  triaxinl  data  to  higher  confining 
pressures  at  a  strain  rate  of  about  JC*^  s’*.  these  data  arc  shown  in  Figure  19. 

5.  msausiQas 

I.  An  approximate  physical  model  for  damage  evolution  in  brittle  solids 
under  compressive  stress  states  has  been  developed.  The  model  is  based  on  the 
growth  of  wing  cracks  from  a  population  of  small,  inclined,  starter  cracks;  and 
the  interaction  between  them.  The  important  variables  of  the  problem  are:  the 
size.  2a,  of  the  initial  inclined  cracks,  and  the  initial  damage 

D  =  ~x  a3Nv 

The  state  of  the  material  is  measured  by  the  current  value  of  the  damage 

D  =  j7c(l  +  aa)  Nv 

where  I  is  the  length  of  the  wing  cracks. 

A  number  of  simplifying  assumptions  are  made  in  the  current  model.  In 
particular  it  is  assumed  that  the  population  of  initial  cracks  all  have  the  same  size. 
One  consequence  of  this  is  that  the  initiation  surface  (that  is  ,  the  combination  of 
stresses  corresponding  to  the  first  extension  of  the  cracks,  and  thus  the  first 
increment  of  new  damage)  is  linear,  described  by: 

<J1  =  Cl  03  -  CT0 
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Similarly,  damage  itself  is  a  state  variable  so  that  surfaces  of  constant  damage  in 
stress  space  arc  also  linear.  The  surfaces  corresponding  to  final  macroscopic 
fracture  is  not  one  of  constant  damage  (as  often  assumed)  .  The  terminal  damage 
itself  depends  on  the  stress  state,  but  it,  too,  is  well  approximated  by  a  linear 
relationship: 


01  aC203-0C 

where  oc  is  the  simple  compressive  strength. 

The  model  has  been  fitted  to  data  for  '*>cr  of  rocks.  The  process  gives 
physical  insight  into  the  damage  accumulation  and  failure  of  these  materials  in 
compression.  In  particular,  the  fitting  process  leads  to  a  value  for  the  coefficient 
of  friction  across  the  crack  faces;  tnc  size  of  the  initial  flaws;  and  the  initial 
damage  Dq  .  The  failure  process  depends  principally  on  these  variables. 

Curvature  of  the  failure  surface  is  shown  to  depend,  at  least  partly,  on  an 
interaction  between  the  brittle  failure  mechanism,  and  plastic  flow.  Rocks  which 
show  clearly  established  plasticity  at  high  pressures  (marbles,  and  rock  salt,  for 
example)  show  a  brittle  regime  at  low  pressures,  a  transitional  regime  at 
intermediate  pressures  (both  depending  strongly  on  pressure),  and  a  regime  of 
plasticity  at  high  pressures  which  is  independent  of  pressure  itself.  It  is 
noteworthy  that  silicate  ror.ks  such  as  granite,  gabbro,  dunite  and  aplite  show  a 
similar  behavior,  with  the  transition  to  plasticity  dominating  failure  at  confining 
pressures  of  general  order  E/30.  This  transition,  at  first  sight  a  surprising  one,  is 
nonetheless  to  be  expected  at  such  stress  levels  which  arc  roughly  the  theoretical 
shear  strength  of  the  minerals  within  the  rock. 

Several  other  noteworthy  conclusions  emerge.  One  is  that,  in  rocks  which 
arc  almost  fully  dense,  the  initial  flaw  size  is  roughly  equal  to  the  grain  size  of  the 
rock  itself.  But  the  initial  damage  level,  Do  ,  varies  widely.  In  law  porosity 
crystalline  silicates  such  as  granite,  this  level  is  low  (typically  3%);  but  in 
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intrinsically-plastic  materials  like  calcitc  ant1  rock  salt,  the  iniliat  damage  level  is 
high  (of  order  15%)  perhaps  because  the  flaws  from  which  wing  cracks  grow  arc 
slip  bands  within  suitably  oriented  grains  rather  than  cracks. 

3.  Data  for  rocks  which  are  almost  fully  dense  arc  well  fitted  by  the  model. 
We  find,  too,  that  porous  rocks  (limestone  and  sandstone,  both  with  roughly  15% 
porosity)  arc  also  well  fitted.  This  suggests  that  an  analogous  theoretical 
development  may  be  possible  for  porosity  induced  cracking  too. 
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Figure  1. 
Figure  2. 

Figure  3. 

Figure  4. 

Figure  5. 

Figure  6. 

Figure  7. 
Figure  8. 
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Compressive  failure  of  a  bridle  solid  containing  a  distribution  of 
flaws. 

Cracks  can  initiate  at  inclined  flaws  and  at  holes.  In  both  eases  there 
arc  two  contributions  to  Ki ,  the  opening  stress  intensity  at  the  tip  of 
the  growing  wing  cracks.  One  is  caused  by  the  stress  concentrations 
at  the  flaw;  the  other  Is  due  to  the  closing  effect  of  03. 

Data  for  crack  initiation  in  Westerly  granite.  Crack  initiation  data 
for  several  other  rocks  arc  analyzed  in  section  4.  In  alt  eases  the 
data  arc  well  fitted  by  equation  (1)  with  a  coefficient  of  friction 
between  O.SS  and  0.65. 

A  population  of  growing  cracks.  We  first  analyze  the  growth  of  an 
isolated  crack  (shown  above)  and  then  include  the  crack-crack 
interaction  (illustrated  on  the  right). 

A  comparison  of  the  approximate  equation  (8b)  with  the  numerical 
calculations  of  Ncmat-Nasscr  and  Horii  (1982).  The  approximation  is 
adequate  for  the  present  purposes. 

The  dependence  of  axial  stress  on  damage  as  predicted  by  equations 

(16)  and  (17).  The  peak  stress  is  marked.  We  take  this  as  the  failure 

stress. 

The  peak  value  of  or]  plotted  against  93  to  show  the  failure  surface. 

In  the  present  formulation,  damage  is  a  state  variable. 

Wing  cracks  growing  from  an  initial,  constrained,  penny-shaped 
flaw.  The  geometry  is  more  complicated  than  in  the  2-dimensional 
case  but  the  same  method  can  be  used  to  give  an  approximate 
solution  for  Kj  at  the  tips  of  the  wing  cracks. 

The  dependence  of  axial  stress  on  damage  as  predicted  by  equations 

(27)  and  (28).  The  peak  stress  is  marked.  We  take  this  as  the  failure 

stress. 

The  peak  value  of  91  plotted  against  93  to  show  the  failure  surface. 

In  this  formulation,  damage  is  a  state  variable. 

Comparison  between  experimental  and  theoretical  failure  surfaces 
for  granite.  Data  and  theory  for  microfracturc  initiation  and 
surfaces  of  constant  damage  are  also  compared.  The  yield  surface 
(cqn.  30)  is  also  plotted  as  the  heavy  broken  lines.  Part  (a)  is  limited 
to  low  and  intermediate  confining  stress;  part  (b)  extends  to  the 
largest  measured  confining  stress. 

Comparison  between  experimental  and  theoretical  failure  surfaces 
for  aplitc.  The  heavy  solid  line  is  the  theoretical  failure  surface. 

The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  arc  surfaces  of 
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constant  damage.  The  heavy  broken  line  is  the  yield  surface  (cqn. 
30).  Part  (a)  is  limited  to  low  and  intermediate  confining  stress;  part 
(b)  extends  to  the  largest  measured  confining  stress. 

Figure  13.  Comparison  between  experimental  and  theoretical  failure  surfaces 
for  dunite.  The  heavy  solid  line  is  the  theoretical  failure  surface. 

The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  arc  surfaces  of 
constant  damage.  The  heavy  broken  line  is  the  yield  surface  (cqn. 
30). 

Figure  14.  Comparison  between  experimental  and  theoretical  failure  surfaces 

for  cclogitc.  Tire  heavy  solid  line  is  the  theoretical  failure  surface. 
The  light  solid  line  is  the  surface  for  the  initiation  of 
microrracturing  while  the  light  broken  lines  arc  surfaces  of 

constant  damage.  The  heavy  broken  line  is  the  yield  surface  (cqn. 
30). 

Figure  IS.  Comparison  between  experimental  and  theoretical  failure  surfaces 

for  gabbro.  The  heavy  solid  line  is  the  theoretical  failure  surface. 
The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  are  surfaces  of 

constant  damage.  The  heavy  broken  line  is  the  yield  surface  (cqn. 
30). 

Figure  16.  Comparison  between  experimental  and  theoretical  failure  surfaces 

for  sandstone.  The  heavy  solid  line  is  the  theoretical  failure  surface. 
The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  arc  surfaces  of 

constant  damage.  The  heavy  broken  line  is  the  yield  surface  (cqn. 
30). 

Figure  17.  Comparison  between  experimental  and  theoretical  failure  surfaces 

for  limestone.  The  heavy  solid  line  is  the  theoretical  failure  surface. 

The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  arc  surfaces  of 

constant  damage.  The  heavy  broken  line  is  the  yield  surface  (cqn. 
30). 

Figure  18.  Comparison  between  experimental  and  theoretical  failure  surfaces 

for  marble.  The  heavy  solid  line  is  the  theoretical  failure  surface. 
The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  arc  surfaces  of 

constant  damage.  The  heavy  broken  line  is  the  yield  surface  (cqn. 

30).  Parts  (a)  and  (b)  are  limited  to  low  and  intermediate  confining 
stress;  part  (c)  extends  to  the  largest  measured  confining  stress. 

Figure  19.  Comparison  between  experimental  and  theoretical  failure  surfaces 

for  rock  salt.  The  heavy  solid  line  is  the  theoretical  failure  surface. 
The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  arc  surfaces  of 

constant  damage.  The  heavy  broken  line  is  the  yield  surface  (cqn. 

30).  Part  (a)  is  limited  to  low  and  intermediate  confining  stress;  part 
(b)  extends  to  the  largest  measured  confining  stress. 
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Figure  20.  Comparison  belween  experimental  and  theoretical  failure  surfaces 
for  marble  at  low  and  intermediate  confining  stress.  The  heavy  solid 
line  is  the  theoretical  failure  surface.  The  light  solid  line  is  the 
surface  for  the  initiation  of  microfracturing  while  the  light  broken 

lines  arc  surfaces  of  constant  damage.  The  heavy  broken  line  is  the 
yield  surface  (cqn.  30). 

Figure  21.  Comparison  between  experimental  and  theoretical  failure  surfaces 

for  marble  at  low  and  intermediate  confining  stress.  The  heavy  solid 
line  is  the  theoretical  failure  surface.  The  light  solid  line  is  the 
surface  for  the  initiation  of  microfracturing  while  the  light  broken 

lines  arc  surfaces  of  constant  damage.  The  heavy  broken  line  is  the 
yield  surface  (cqn.  30). 

Figure  22.  Comparison  belween  experimental  and  theoretical  failure  surfaces 

for  marble  extended  to  the  highest  measured  confining  stress.  The 
heavy  solid  line  is  the  theoretical  failure  surface.  The  light  solid 
line  is  the  surface  for  the  initiation  of  microfracturing  while  the 
light  broken  lines  are  surfaces  of  constant  damage.  The  heavy 
broken  line  is  the  yield  surface  (cqn.  30). 

Figure  23.  Comparison  between  experimental  and  theoretical  failure  surfaces 

for  rock  salt  at  low  and  intermediate  confining  stress.  The  heavy 
solid  line  is  the  theoretical  failure  surface.  The  light  solid  line  is 

the  surface  for  the  initiation  of  microfracturing  while  the  light 

broken  lines  arc  surfaces  of  constant  damage.  The  heavy  broken 

line  is  the  yield  surface  (eqn.  30). 

Figure  24.  Comparison  between  experimental  and  theoretical  failure  surfaces 

for  rock  salt  extended  to  the  largest  confining  stress.  The  heavy 
solid  line  is  the  theoretical  failure  surface.  The  light  solid  line  is 

the  surface  for  the  initiation  of  microfraciuring  while  the  light 

broken  lines  are  surfaces  of  constant  damage.  The  heavy  broken 

line  is  the  yield  surface  (cqn.  30). 
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Compressive  failure  of  a  brittle  solid  containing  a  distribution  of  flaws. 


J  PI  I  I  PI  I 


Figure  2.  Cracks  can  initiate  at  inclined  (taws  and  at  holes.  In  both  cases  there  are 
two  contributions  to  K| .  the  opening  stress  intensity  at  the  tip  of  the 
growing  wing  cracks.  One  is  caused  by  the  stress  concentrations  at  the 
(law;  the  other  is  due  to  the  closing  effect  of  03. 
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Figure  3.  Data  for  crack  initiation  in  Westerly  granite.  Crack  initiation  data  for 

several  other  rocks  are  analyzed  In  section  4.  In  all  cases  the  data  are  well 
fitted  by  equation  (1)  with  a  coefficient  of  friction  between  0.55  and  0.65. 
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A  population  of  growing  cracks.  We  first  analyze  the  growth  of  an  isolated 
crack  (shown  above)  and  then  include  the  crack-crack  interaction 
(illustrated  on  the  right). 
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Figure  6.  The  dependence  of  axial  strew  on  damage  as  predicted  by  equations  (16) 
and  (17).  The  peak  stress  Is  marked.  We  take  this  as  the  failure  stress. 
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Figure  7. 


The  peak  value  of  01  plotted  against  03  to  show  the  failure  surface.  In  the  ' 
present  formulation,  damage  is  a  state  variable. 
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Figure  8.  Wing  cracks  growing  from  an  initial,  constrained,  penny-shaped  flaw.  The 
geometry  is  more  complicated  than  in  the  2-dimensionai  case  but  the  same 
method  can  be  used  to  give  an  approximate  solution  for  K|  at  the  tips  of  the 
wing  cracks. 
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Figure  9.  The  dependence  of  axial  stress  on  damage  as  predicted  by  equations  (27) 
and  (28).  The  peak  stress  is  marked.  We  take  this  as  the  failure  stress. 
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Figure  10.  The  peak  value  of  ai  plotted  against  03  to  show  the  failure  surface.  In  this 
formulation,  damage  is  a  state  variable. 
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AXIAL  STRESS,  o',  (MRq) 


-800  -600  -400  -200  0 


Figure  11.  Comparison  between  experimental  and  theoretical  failure  surfaces 
for  granite  at  low  and  intermediate  confining  stress.  Data  and 
theory  for  microfracture  initiation  and  surfaces  of  constant  damage 
are  also  compared.  The  yield  surface  (eqn.  30)  is  also  plotted  as  the 
heavy  broken  lines. 
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RADIAL  STRESS,  cr3  (MPa) 


AXIAL  STRESS  cr,  (MPa) 


Figure  12.  Comparison  between  experimental  and  theoretical  failure  surfaces 
for  granite  extended  to  the  largest  measured  confining  stress.  Data 
and  theory  for  microfracture  initiation  and  surfaces  of  constant 
damage  are  also  compared.  The  yield  surface  (eqn.  30)  is  also  plotted 
as  the  heavy  broken  lines. 
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AXIAL  STRESS,  c rt  (MPa) 


Figure  13.  Comparison  between  experimental  and  theoretical  failure  surfaces 
for  aplite  at  low  and  intermediate  confining  stress.  The  heavy  solid 
line  is  the  theoretical  failure  surface.  The  light  solid  line  is  the 
surface  for  the  initiation  of  raicrofracturing  while  the  light  broken 
lines  are  surfaces  of  constant  damage.  The  heavy  broken  line  is  the 
yield  surface  (eqn.  30). 
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AXIAL  STRESS,  <r,  (MPa) 

-4000  -2000  0 


Figure  14.  Compirison  between  experimenul  and  theoretical  failure  surfaces 
for  aplite  extended  to  the  largest  measured  confining  stress.  The 
heavy  solid  line  is  the  theoretical  failure  surface.  The  light  solid 
line  is  the  surface  for  the  initiation  of  microfracturing  while  the 
light  broken  lines  are  surfaces  of  constant  damage.  The  heavy 
broken  line  is  the  yield  surface  (eqn.  30). 
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RADIAL  STRESS,  cr,  (MPa) 


AXIAL  STRESS,  o;  (MPa) 


NORMALIZED  AXIAL  STRESS,  q  /KxlO-3 


Figure  15.  Comparison  between  experimental  and  theoretical  failure  surfaces 
for  dunite.  The  heavy  solid  line  is  the  theoretical  failure  surface. 
The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  are  surfaces  of 
constant  damage.  The  heavy  broken  line  is  the  yield  surface  (eqn. 
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RADIAL  STRESS,  cr  (MPa) 
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AXIAL  STRESS,  cr,  (MPa) 
-6000  -3000 


NORMALIZED  AXIAL  STRESS,  q/KxlO"3 


Fiwre  16.  Comparison  between  experimental  and  theoretical  failure  surfaces 
for  eclogite.  The  heavy  solid  line  is  the  theoretical  failure  surface. 
The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  are  surfaces  of 
constant  damage.  The  heavy  broken  line  is  the  yield  surface  (eqn. 
30). 
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Figure  17.  Comparison  between  experimental  and  theoretical  failure  surJ**;*s 
for  gabbro.  The  heavy  solid  line  is  the  theoretical  failure  surface. 
The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  are  surfaces  of 
constant  damage.  The  heavy  broken  line  is  the  yield  surface  (eqn. 

30). 
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Fiffure  19  Comparison  between  experimental  and  theoretical  Mure  iurfi^5s 
F““re  W-  fer  limcILc.  The  be.vy  .olid  lino  i.  <he  theoretical  f..lure  surface. 
The  light  solid  line  is  the  surface  for  the  initiation  of 
raicrcfncturing  while  the  light  broken  lines  are  Jurfacesof 
constant  damage.  The  heavy  broken  line  is  the  yieio  surface  (eqn. 

30). 
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Figure  20.  Comparison  between  experimental  and  theoretical  fail 

for  marble  at  low  and  intermediate  confining  *tress;  yhc  I1  X 

line  is  the  theorerical  failure  surface.  The  light  solid  I line  is  th ic t 
surface  for  the  initiation  of  raicrofractunng  white  the  light  broken 
lines  are  surfaces  of  constant  damage.  The  heavy  broken  line  is  the 

yield  surface  (eqn.  30). 
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Figure  21.  Comparison  between  experimental  and  theoretical  tailure  surfaces 
for  marble  at  low  and  intermediate  confining  stress.  The  heavy  solid 
line  is  the  theoretical  failure  surface.  The  light  solid  line  is  the 
surface  for  the  initiation  of  microfracturing  while  the  light  broken 
lines  are  surfaces  of  constant  damage.  The  heavy  broken  line  is  the 
yield  surface  (eqn.  30). 
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Figure  22.  Comparison  between  experimental  and  theoretical  failure  surface; 

for  marble  extended  to  the  highest  measured  confining  stress.  The 
heavy  solid  line  is  the  theoretical  failure  surface.  The  light  solid 
ljne  is  the  surface  for  the  initiation  of  microfracturing  while  the 
light  broken  lines  are  surfaces  of  constant  damage,  TCie  heavy 
broken  line  is  the  yield  surface  (eqn.  30). 
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Figure  23.  Comparison  between  experimental  and  theoretical  failure  surfaces 
for  rock  salt  at  low  and  intermediate  confining  stress.  The  heavy 
solid  line  is  the  theoretical  failure  surface.  The  light  solid  line  is 
the  surface  for  the  initiation  of  microfracturing  while  the  light 
broken  lines  are  surfaces  sf  constant  damage.  The  heavy  broken 
line  is  the  yield  surface  (eqn.  30). 
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Figure  24.  Comparison  between  experimental  and  theoretical  failure  surfaces 
for  rock  salt  extended  to  the  largest  confining  stress.  The  heavy 
solid  line  is  the  theoretical  failure  surface.  The  light  solid  line  is 
the  surface  for  the  initiation  of  microfracturing  while  the  light 
broken  lines  are  surfaces  of  constant  damage.  The  heavy  broken 
line  is  the  yield  surface  (eqn.  30). 
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